# Full text of "Radial multiresolution in dimension three"

## See other formats

Radial multiresolution in dimension three I Holger Rauhut and Margit Rosier O ■ O ■ (N : _ Abstract ^ . We present a construction of a wavelet-type orthonormal basis for the space of radial ' L^-functions in R'^ via the concept of a radial multiresolution analysis. The elements . of the basis are obtained from a single radial wavelet by usual dilations and generalized (N ■ translations. Hereby the generalized translation reveals the group convolution of radial functions in M'^. We provide a simple way to construct a radial scaling function and a ' radial wavelet from an even classical scaling function on M. Furthermore, decomposition I ■ and reconstruction algorithms are formulated. ! 2000 AMS subject classification: 42C40, 43A62 ^ ' Keywords: wavelets, multiresolution analysis, radial functions, generalized translation, Bessel- G ■ Kingman hypergroup ^ ■ 1 Introduction (N : I Standard approaches in multivariate wavelet analysis are based on the construction of mul- ■ tiresolution analyses and wavelet bases from aflfine transformations of a finite set of basis , functions, called multi- wavelets (see e.g. |13| . [21). The translations are taken from a lattice I subgroup r of {W^, +) and the dilations are given by the integer powers of an expansive ' matrix which leaves T invariant. The number of multi- wavelets needed to obtain a full basis . of L^(M'^) depends on the determinant of the dilation matrix and is in general larger than g ; 1 for d > 1. ■ However, if one restricts to the analysis of radially symmetric functions only, it suggests K»" ■ itself to exploit this symmetry in the construction of corresponding wavelet transforms in ^ I order to reduce the high amount of computational effort. A purely radial setting would ^ ' for example naturally occur when separating variables in polar coordinates, and treating ■ ■ the spherical and radial parts separately. There is a broad literature dealing with wavelet analysis and multiresolution on spheres, see e.g. PP, [HI and the references therein. In the radial case, it is not difficult to establish a continuous wavelet analysis based on the convolution structure of radial functions or measures instead of the usual translation in W^. Radial convolution structures are special cases (for half- integer indices) of so-called Bessel-Kingman hypergroups (see jl], J2]), and a continuous wavelet analysis can in fact be developed in this general setting, see e.g. ^H], and ^5]. Essentially the same concept is underlying the approach of Epperson and Frazier where radial wavelet expansions in are constructed which are based on sampling lattices with the spatial discretization determined by the positive zeros of related Bessel functions of the first kind. The spatial lattice is equidistant only in the special cases d = 1 and d = 3. This can be seen as an obstruction against a multiscale approach to radial wavelets in arbitrary dimensions. 1 As to the authors knowledge, radial multiresolution analyses have in fact not been con- sidered up to now, and there seems to be no general rigorous approach available for the construction of orthogonal radial wavelet bases in arbitrary dimension. This problem is closely related with the question to find a Poisson summation formula compatible with the Bessel-Kingman translation, which still remains open. We mention that the construction of j7j relies on the requirement that the involved radial wavelets are band-limited, i.e. their Fourier transforms have compact support. In the present paper, we construct radial multiscale analyses and orthogonal radial wavelet bases in M^. In dimension 3, the algebraic structure of the radial convolution allows to carry out the constructions along the same lines as in the well-known Euclidean setting. Hereby, the equidistance of the zeros of the corresponding Bessel function (which is simply a sinc-function) is of decisive importance. In order to motivate our approach and to make the intrinsic problems towards an extension to arbitrary dimensions more visible, we start the paper with a short account on the continuous wavelet transform for Bessel-Kingman hypergroups, and also present some material about Bessel frames. This in particular comprises the radial case in arbitrary dimension. The continuous transform is, up to normalization, the same as in ^H], while the Bessel frames are constructed in the spirit of the radial wavelet bases in jjj. In particular, in the spatial discretization the zeros of associated Bessel functions occur in a natural way, and the wavelets are band-limited. The radial analysis in is then also carried out in the setting of the corresponding Bessel-Kingman hypergroup H on [0, oo). Our concept of a radial multiresolution analysis (MRA) in is in fact that of a MRA for the L^-space L'^{H) of this hypergroup. The scale spaces {Vj)j^i C L^{H) are obtained by dyadic dilations from Vq, which in turn is spanned by equidistant hypergroup translates of a fixed "radial" scaling function (p G L^(H). It is characterized by a two-scale relation, but in contrast to classical MRAs, (j) itself is not contained in Vq, and the scale spaces are not shift-invariant (w.r.t. the hypergroup translation). Particular emphasis is put on the construction of orthogonal MRAs. Here periodicity arguments similar to the classical case are needed which would not be available in arbitrary dimensions. From a given orthogonal MRA we then derive an orthogonal wavelet basis for the underlying hypergroup. By construction, this "radial" basis has a direct interpretation as an orthogonal wavelet basis for the subspace of radial functions in L-^(R'^). We also provide a concise characterization of radial scaling functions in terms of even classical scaling functions on R. It in particular implies that in contrast to the classical case, there do not exist any real-valued orthonormal radial scaling functions with compact support. The paper is organized as follows: In Section 2, we recall basic facts from the analysis of radial functions in R*^, explain the corresponding radial hypergroup convolution struc- ture, and extend the setting to Bessel-Kingman hypergroups of arbitrary index. Section 3 contains a short account on the continuous wavelet transform based on the Bessel-Kingman translation, as well as the construction of Bessel frames. In Section 4, radial multiresolution analyses in R^ and their scaling functions are introduced and discussed, while Section 5 is devoted to the construction of orthogonal radial wavelet bases. The connection between radial scaling functions in R^ and classical scaling functions on R is established in Section 6. Finally, in Section 7 decomposition and reconstruction algorithms are discussed. 2 2 Radial analysis and Bessel- Kingman hypergroups Suppose F e L2(M'^) is radial, i.e. F{Ax) = F{x) a.e. for all A e SO{d). Then there is a unique / G L^{M.^,lO(1/2-i) such that F{x) = where | . | denotes the Euclidean norm on and for a > 0, the measure Ua on = [0, oo) is defined by dujair) = (2"r(a + l))~V2°^+^dr. Its normalization implies that ||-F||2 = H/lb, t^^/a-i' '^herc || . ||2 is taken with respect to the normalized Lebesgue measure {2Tx)~'^/'^dx on M"^. On L^(M_|_,a;a) the Hankel transform of index a is defined by /•oo r(A) = / U\r)f{r)dw^{r) Jo with the normalized Bessel function n=0 ^ ' There is a Plancherel Theorem for the Hankel transform, which states that / i— /" es- tablishes a self-inverse, isometric isomorphism of -L^(]R+, Wq,). If F G L^(R'^) is radial with F{x) = /(|a;|), then a short calculation shows that its Plancherel transform is again radial with J^{F){^) = /^'^/^"^Hl^D- This is due to the fact that Jd/2-i(kl)= / e-(^'f>da(0, (2.1) where da denotes the spherical surface measure normalized according to da{S'^~^) = 1. In contrast, the usual group translates x i-^ F{x + y) , y <E M.'^' , will no longer be radial (apart from trivial cases). However, we observe that the spherical means ■F{x) := / F{x + rS,)da{S,), re of F are again radial. Moreover, ||MrF||2 < \\F\\2- Thus Mr induces a norm-decreasing linear mapping Tr : L\R+,ua/2-i) ^ L\R+,ua/2-i), Trfi\x\) := MrF{x), where / and F are related as above. Put a = d/2 — 1. Then a short calculation in polar coordinates gives Trfis) = Ca r /(Vr2 + s2-2rscos^) sin^'^ipdip with = ^/^"^w^ - (2-2) Jo 1 (a-l- 2ji (2) This defines a norm-decreasing generalized translation on L'^(R^,uOa) not only for a = d/2—1, but also for general a > —1/2. Having harmonic analysis in mind, we are thus lead 3 to introduce a corresponding measure algebra on M+r For r,s £ we define a probability measure 6r *a on M+ by Sr *a 6sU') ■■= Ca C f {^f + - 2rs COB ^) sin2°(^d(^, / G C,(M+). (2.3) (Cc(K_|-) denotes the space of continuous, compactly supported functions on M_|_). The con- volution (|2.3j) of point measures extends uniquely to a bilinear, commutative, associative and weakly continuous convolution on the space M{,(M_(_) of regular bounded Borel measures on M_|_. It is probability-preserving and makes -/Vfb(M_|_) a commutative Banach-*-algebra with respect to total variation norm, with neutral element 5q and the mapping fx >-^Jl as involu- tion. The pair (M+,*q,) is called the Bessel-Kingman hypergroup of index a. We write instead of M+ when putting emphasis on the specific convolution structure. Generally speak- ing, a hypergroup is a locally compact Hausdorff space together with a weakly continuous and probability preserving convolution of regular bounded Borel measures generalizing the measure algebra of a locally compact group; it also has a unit and an involution substituting the group inverse. In particular, every locally compact group is also a hypergroup. There is a well-established harmonic analysis for commutative hypergroups, which in the special case of reflects the harmonic analysis of radial functions (and measures) in M*^. We refer the reader to or |1] for a general background on hypergroups, including ^] for the Bessel-Kingman case. Let us mention only some aspects which are of importance in our context: The measure uJa is a Haar measure for H^, i.e. it satisfies / Tsfduja= / fdu^ for ah / G C7,(M+). Jo Jo Up to a constant factor, uJa is the unique positive Radon measure on M_|_ with this property. Moreover, / {Tsf)gduj^= / fiTsg)duJa for s G R+ (2.4) Jo Jo whenever both integrals exist. The Bessel functions satisfy the product formula 5r *a Ss{ja) = ja{r)ja{s) for all r, s G M+; see |171 11.4]. For half- integers a, this is easily deduced from 1)2.1(1 . In fact, the functions f jai^i"), A G M+ are exactly those which are bounded and multiplicative with respect to *a- They constitute the so-called dual space of the hypergroup Ha- In this way, the Hankel transform / i— > /° on L^(M+,a;Q) can be interpreted as a Plancherel transform for Ha- For abbreviation, we put U'{Ha) := LP(M+,Wa) and we denote by (., .) and || . ||2 the scalar product and norm in L'^iHa)- It follows easily from ((TH) that for / G L'^{Ha), xTfiX) = ja{Xr)riX)- (2.5) 3 Continuous wavelet transform and frames for Bessel-Kingman hypergroups In order to put the multiresolution approach in the following sections into a suitable frame- work, we continue with a short account on the continuous wavelet transform and wavelet frames for Bessel-Kingman hypergroups. This in particular includes a continuous radial wavelet transform and radial wavelet frames in arbitrary dimensions. 4 3.1 The continuous wavelet transform The following construction is essentially the same as in |16j , only with a different normaliza- tion of the dilation operators and the resulting wavelet transform. (In contrast to we choose dilations to be unitary, see below). We shall therefore be brief in our presentation, and refer the reader to [IBj Sect. 6. Ill] for further details. Let B{L'^{Ha)) denote the space of continuous linear operators on L'^{Ha)- Besides the translation operators Tr £ B{L?'{Ha)) introduced in Section 2, we consider the dilations Daf{r) «>0 which are obviously unitary in L'^(Ha)- Notice also that (AJ)" = DyJ^. (3.1) We define TT-.H^x (0, oo) B{L\Ha)), 7r(r, a) := Tr-Da- It is easily checked that vr is continuous with respect to the weak operator topology on B{L'^{Ha)), c.f. Prop. 6.III.7]. Definition 3.1. A function g £ L^{Ha) is called admissible, if For abbreviation, put diJa{r,a) := ^2a+3 daduja{r). The following is a reformulation of |161 Thm 6.III.1] in terms of our notation: Theorem 3.2. (Plancherel Theorem) If g G L'^{Ha) is admissible, then for all f G L'^{Ha), [ K/,^(r,a)5)|'a!5,(r,a) = • ||/||2 . JHaX{0,oo) Polarization further implies for admissible gi, 52 and arbitrary /i, /2 S L^{Ha) the orthogonality relation dX A {fi,7r{r,a)gi){f2,7r{r,a)g2)duJa{r,a) = (/i,/2)- / 5f(A)?2°(-^) 'Ha X (0,00) Jo Definition 3.3. Let g € L'^{Ha) be admissible. The mapping $g : L^{Ha) — > L^{Ha X (0,oo),5„), ^'g/(r,a) := {f,7r{r,a)g) is called the wavelet transform on Ha with analyzing wavelet g. For a = d/2 — \, it coincides with the continuous wavelet transform on of a radial function f with radial wavelet g, see rr^ . Some types of inversion formulas for this transform can be found in ^ 5 3.2 Bessel frames Let us now turn to possible discretizations. In order to obtain discrete versions of the usual wavelet transform on M, it is standard to use sampling lattices of the type {{nba^, a^), k,n ^ Z} with constants a > 1, 6 > 0. Here the discretization of the translation parameter is in accordance with the related group structure of M and is therefore (in general) not appropriate for radial wavelet analysis. Following Epperson and Frazier [Jj, we propose lattices where the discretization of the translation parameter involves the positive zeros Q < Va,i < '^a,2 < ■ ■ ■ of the Bessel function ja- By a result of McMahon, these are asymptotically distributed according to A standard lattice in Ha x (0, cxd) is given by { {va,nha^ ,a^), A; G Z, n G N}, (a > 1, 6 > 0). The "almost orthogonal" radial wavelet expansions of Epperson and Frazier [TJ are based on this type of sampling lattice (with a = 2). In the following result, the discretization of the dilation parameter is still rather arbitrary. Theorem 3.4. LetQ be a countable subset of {0,oo) and g G L?'{Ha)- Assume that suppg°' is contained in [0, /] for some I > 0, i.e. = a.e. on (/, cxd), and that there exist constants A, 5 > such that -4 < ^ |?"(^A)p < B for almost all A G M+ . For n G Z and q € Q define "wavelets" gn,q G L?'{Ha) by gn,q ■■= ■ Tr^qDq{g) = 7r{rnq,q)g, where 1 -V, and M" 2(i-")/2 J," t Vr(a + 1) I J„+i(z^„,„j| Then the set {gn,q '■ n. G N, g G Q} is a frame for L?'{Ha) with frame bounds A/^"^^ and BZ2"+2. This means that for f e L'^{Ha), Ai'-^'-\\m<EE\(9n,,j)\' <Bi'"^'-\\f\\i qeQneN Proof. The decisive point in the proof is the fact that the normalized Fourier-Bessel functions p^{X) := M„"j„(z/„,„A), nGN form an orthonormal basis of the Hilbert space Xa := ^^([0, 1], cjalfo,!]); see e.g. Erdelyi et al. Using (|TT|) and (^31) . we write ?„",(A) = M^Ur^qX)Dyqg^iX) = p^{^-X)Dyqg"{X). 6 By the Plancherel theorem for the Hankel transform, we obtain where for the last identity, we used that the support of Di/i{'g°') is contained in [0,1]. Parseval's identity for now yields pi poo T.\(9n,gJ)\'= / \Dyi{g")D^/i{f-)fdcOa = / \g''{qX)\' \r{X)\' dua{X). _^,>,T JO Jo Hence, poo , , This implies the assertion. □ The second condition of this theorem is rather implicit. Following e.g. Bernier and Taylor [3, it is possible to obtain sufficient criteria which are easier to check by introducing the concept of "frame generators" . But in order to stay concise, we restrict ourselves to the most interesting special case of a standard lattice as defined above. Here Q = {a*^. A; G Z} with a > 1. Proposition 3.5. Suppose that g G Lp'^Ha) satisfies the following conditions: (i) has compact support which is contained in the open interval (0,oo). (ii) essinf\\g°'{\)\ : A G [a^jO""*"^]} > a for some n G Z and a > 0; (iii) r := ||5°||oo < oo. Then there exists a constant M > such that CT^ < J]] |5°(a''A)|2 < Mr^ for almost all A G M+ . (3.2) fcez Consequently, the set {gna'' '■ G N, A; G Z} is a frame for L^{Ha) with bounds cr^Z^"'''^ and Mr2/2"+2. Proof. We use the arguments of [3], sect. 4 in a simplified form which is adapted to our situation. Put F := [l,a). Then the intervals a^F = [a"',a^~^^), n G form a disjoint cover of (0, oo). As T := supp^" is compact in (0, oo), it is covered by finitely many of the a"F. This implies that M := sup ^{keZ: AG a^'^T} Ag(0,oo) is finite. By (iii), this gives the upper bound in (|3.'2j) . The lower bound follows from (ii) together with the fact that the a"'F cover (0, oo). □ 7 4 Radial Mult iresolut ion Analysis in Radial analysis in corresponds to the Bessel-Kingman hypergroup Ha with a = 1/2. For convenience we shall usually omit the subscript 1/2 and put duj{r) := duJi/2{r) = J -r'^dr, f := f^'^, j{r) := ji/2(r) = ^— . We further write H instead of -ffi/2 a.nd denote by ( . , . ) and || . ||2 the scalar product and norm in L?'{H) = L'^{H,ij), respectively. Notice that the Bessel function j is even on M. Hence it is natural to assume the Hankel transform /•oo fix) = / jiXr)f{r)du;ir) Jo of / G Lp'iH) to be continued to an even function on M as well. We shall always do this throughout the paper. We also mention that by a change of variables, the generalized translation (|2.2|) on H can be written in the simple form Trfis) = ^ f{t)tdt. (4.1) The non-negative zeros of the Bessel function j are given by tk := kiT, k e'N , and the normalized Fourier-Bessel functions Pk{r) := Mkjihr) with = 2^/\''/^k (4.2) form an orthonormal basis of the Hilbert space L^([0, l],u;|[o,i])- This is equivalent to the obvious fact that the functions s^{r) := (2/7r)^/^ rpfc(r) = \/2 sin(/c7rr), fcGN are an orthonormal basis for [0, 1] := -^^^([0, 1] , dr). It will be of importance in the following that the Sk are 2-periodic. Let us come to the definition of a radial multiresolution analysis (MRA) for W^, i.e. for the Bessel-Kingman hypergroup H. It is close to the well-known definition of Mallat Jl] for M. For convenience, we introduce the notation If / G L'^{H) then according to (^3]) . (M,rW/)-(A) = p,(A)/(A) = ilf-^fiX). (4.3) Definition 4.1. (Radial Multiresolution Analysis) A radial MRA /or M'^ is a sequence {Vj}j<^i of closed linear subspaces of L^{H) such that (1) Vj C Vj+i for all j G Z; 8 (2) nr=_oo = {0}; (3) Uj°l-oo is dense in L\H); (4) / G Vj if and only i//(2-) G Vj+i; (5) There exists a function 4> G such that := {MkT^^U -.ken} is a Riesz basis ofVo, i.e. spanB^ is dense in Vq and there exist constants A,B>0 such that < ||X]afcMfcTW(/<||i < B\\a\\l k=l for all a = {ak)km G ^^(N); here \\a\\2 = (Zl^i lafcP) The function (p in (5) is called a scaling function for the MRA {V^ }. We remark explicitly that in contrast to the classical case, (p itself is not contained in Vq and Vq is not shift invariant; in fact, if / G Vq then T^'^^f ^ Vq for all k. This will be shown in Corollary 14.61 below. Our first aim is to determine an orthonormal basis for Vq from its Riesz basis, i.e. a function cj)* G L^{H) such that B^p* consititutes an orthonormal basis for Vq. For (p G Lp'{H) we define oo P^A) := 1^(^ + 2^)1'' n=— oo which is even and 2-periodic on M. Proposition 4.2. Let (j) G L'^{H) and A,B >Q. Then oo ^llalli < II akMkT^^^ \\l < B\\a\\l for all a G /^(N) (4.4) k=l if and only if A<P^{\)< B for almost all A G M. (4.5) Proof. Let a G /^(N) be an arbitrary finite sequence. Define oo 5:=^QfcSfc gL2[0,1]. (4.6) k=l We may regard a as an odd, 2-periodic function on M. By the Plancherel theorem for the Hankel transform and ()4.3() . oo oo oo ^QfcMfcTW(/<||2 = ll^afcpfc^lli = / |^afcSfc(A)|2|^(A)|2dA 7 1 / 1 7 1 k=l k=l k=l |5(A)|2|(/>(A)|2dA = - I J5(A)|2p^(A)dA = / |5(A)|^P4A)dA. The Sfc forming an orthonormal basis of L^[0, 1], we have ||a!||2 = ||ct||£^2[o,i]- As the finite sequences form a dense subspace of /^(N), this implies the assertion. □ 9 With A = B = 1 we immediately obtain Corollary 4.3. For (j) G L^{H) the following statements are equivalent: (1) The set = {MkT^^^ : G N} is orthonormal in L'^{H); (2) P^ = l a.e. For 4> G L'^{H), put := span_B<^ , the closure beine; taken in L^{H). The set B^ is a Riesz basis of if and only if there exist constants A,B > such that the equivalent conditions of Proposition are satisfied. This will be a standard requirement in the sequel, and we therefore introduce a separate notation: Definition 4.4. A function <p £ satisfies condition (RB) if B^ is a Riesz basis of As before, we shall often consider functions from L^[0, 1] as odd, 2-periodic functions on M. We therefore define S := {a : M — > C | a|[o_i] G ^^[0, 1], a{—x) = — a(x), a{x + 2) = a{x) for almost all x}. is a Hilbert space with norm || . ||l2[o_i]. Lemma 4.5. Let cp G satisfy (RB). Then for f G Lp'{H) we have the equivalence feV^^ /(A) = ^0(A) with (3eS. The function / G corresponding to (5 = Yl'kLi '^kSk £ S with {ak)keN G /^(N) is given by Proof. By (|4.3() we have ^^(A) = {^f\MkT(''UrW. The translates M^T^^^cj) form a Riesz basis of V^, and hence the (Mf^T^^^cf))^ are a Riesz basis of V(j,. As {sk)km is an orthonormal basis of S, this implies the assertion. □ Lemma 14.51 is of particular interest when (p is the scaling function of a MRA {Vj}. Then Vq = V^, and we easily deduce the previously mentioned lack of shift-invariance: Corollary 4.6. Let {Vj}j(zz be a radial MRA. Then f e Vq implies that T'^'^^f = TtJ i Vq for all G N. Similarly, f & Vj implies that T2-jt^ f ^ Vj for all A; G N. Proof. After rescaling it is enoug h to consider Vq. Recall that Mk{T^''^f )^ = pkf. But if f3 £ S, then pkP ^ S for all k, because periodicity is lost. The characterization of = Vq according to the previous Lemma thus shows that for / G Vq, T^^'^f ^ Vq. □ In the situation of the Lemma, we can easily determine an orthonormal basis of by renormalization: 10 Theorem 4.7. (Orthogonalization) Suppose (j) G L?'{H) satisfies condition (RB). Define (jf G LP'{H) by its Hankel transform Then B^* = {M^T^^^^cp* : G N} forms an orthonormal basis ofV^ = V^* . If (/> is a scaling function of a MRA {V^ }, then = V^* = Vq, and we call (j)* an orthonormal scaling function for {Vj}. Proof. By definition of 0* we have P,^* = 1 a.e. and hence i?,^* is orthonormal according to Corollary 14.31 It remains to prove that V^* = V^. For this, we have to verify that MfcT^'^V* G V4, and MkT^^^ G V^. for all /c G N. Employing Lemma ESI relation and finally the relation Sfc(r) = (7r/2)-'^/^ rpfc(r), one obtains that the above conditions are equivalent to G S, Sky^ G S for all A; G N. But these conditions are obviously satisfied by our assumption on . □ Let us return to our definition of a radial MRA for M'^. Suppose we start with a function (j) G L^{H) satisfying condition (RB) with Riesz constants A,B>0. Define corresponding scale spaces {Vj}ji^z by where the dilation operator Da G B{L^{H)) is defined as in Section 3. Then in particular, the Vj satisfy axiom (4) of Definition 14.11 Put further 0,-fc(r) := D2-, {MkT^''U){r) = Mk{T^^^ (t>){2^ r) j£Z,kGn. (4.8) Then {(pj^k, = {4>o,k-, 4'o,i) for all j, k, I. Thus the {(/'j,^, A: G N} form a Riesz basis of Vj, with the same Riesz constants A, B as for j = 0. In particular, Vj = span{4>j^k,k G N}. Moreover, if B^ = {(po^k '■ A; G N} is an orthonormal basis for Vq, then {(pj^k, A; G N} is an orthonormal basis of Vj. Recall now axiom (1) of Def. 14.11 which requires that the Vj are nested. Similar as in the classical case, this condition can be reformulated in terms of a two-scale relation for (p: Proposition 4.8. For cp with (RB) and {Vj}j,^z as above, the following statements are equivalent: (1) Vj C Vj+i for all j G Z (2) F-i C (3) There exists a function 7 G such that sin(27rA)^(2A) = 7(A)^(A). (4.9) 11 In this case, the coefficients {hk)keN £ ^^(N) in the two-scale relation oo <^-i,i = ^hk(^o,k (4.10) k=l are the coefficients in the Fourier sine series 0/7 G S: 1 °° 00 i 7 2 k=l Proof. Rescaling by the factor 2^ shows that (1) and (2) are equivalent. For (2), we need at least (/>-i,i G Vq. According to Lemma 14.51 this is equivalent to the existence of a function (3 ^ S such that for almost all A, \/8Mij(2vrA)^(2A) = ^i(A) = ^0(A). (4.11) Moreover, if has the expansion (|4.1flj) . then (3 = {tt /2Y^^ Yl^=i ^kSk- In turn, equation ()4.11() is equivalent to relation 1)4. 9() with 7 = (87r)-i/4/?. This gives the stated connection between 7 and It remains to show that € Vq (or equivalently, relation (|4.9|) ) already implies that 0_i^fc G Vq for all A; G N. Similar as above, the latter is equivalent to sin(2A;7rA)^(2A) = A;7fc(A)^(A) (4.12) with 7^, G S. The relation between 7^ and 0-i fc is now given by 7fc = -^^^f^^k, 4>-i,k = ^h\^''4>o,i- 1=1 1=1 Comparison of (|4.9|) with (|4.12l) yields 7fe(A) = iW^^^^^^, = 7(A)^fc-i(cos2vrA) (4.13) k sm(27rAj where sin(A; + l)t Uk{x) = , . . , x = cost (k + 1) sm t denotes the fe-th Chebychev polynomial of the second kind, normalized such that Uk{l) = 1. Thus given 7 G 5", we define 7fc(A) := 7(A)[/fc_i(cos27rA). As t/fc-i is bounded on [—1, 1], 7^ is contained in S as well, and hence (p-i^k ^ ^o- D Let us now consider the remaining axioms (2) and (3) of a radial multiresolution analysis. Theorem 4.9. Let (p G Lp'{H) satisfy condition (RB) and assume that the scale spaces Vj = span{(pj^k : fc G N}, i G Z satisfy V-i C Vq. Suppose further that \(f>\ is continuous in 0. Then {Vj}j^z is a radial MRA if and only if 0(0) 7^ 0. Moreover, <j) is an orthonormal scaling function if and only if \m\ = 1- 12 We remark that continuity of </) in (even on is for example guaranteed if (/> G L'^iH)nL^H). Proof. We have to check axioms (2) and (3). This may be done by shght modifications of standard arguments in the affine case. The condition on (p in will be needed only for (3). We define an orthonormal scaling function (f)* according to Theorem 14.71 The orthogonal projection Pj of L'^{H) onto Vj is then given by klVj,k where the (jf- ^ are defined as in (|4.8j) . For (2), we need to show that limj^_oo H-Pj/lb = for all / G LP'{H). Since functions with compact support are dense in L'^{H), we may assume that supp/ is contained in a compact interval [0, ii]. Parseval's equation then implies k=l k=l fR °° (■2m °° <\\f\\l- J2\^lk{r)\'du;{r) = \\fg- ^\cPl,{r)\'dco{r). ■^0 k=l fc=l Using the explicit formula (|4.1j) for the hypergroup translation in H we further deduce -(1/4 fkn+r 2 _ _ (l)*{t)tdt k=l ' k=l fc=l' J\kn-r\ 00 00 00 /r) xl/4 fkn+r Now assume that j is sufficiently small so that 2^R < tt/2. Then for r E [0,2-^i?], the integration domains [kn — r, /cvr + r] do not overlap and we obtain El'^S..(OP < 7E/ W{t)\H'dt<'' * l|2 — IIV^ 2 r k=l k=l' with suitable constants C, C" > independent of j. Hence for j sufficiently small, \Pjf\\l < C" / -dcoir), Jo which tends to as j ^ — oo. This proves (2). As to (3), suppose first that 0(0) / and let /i G ( U^_^ Vj)^ , i.e. Pjh = for all j £ Z. We claim that h = 0. Indeed, for e > there exists a function / G L?'{H) such that the support of its Hankel transform / is compact and ||/ — /i||2 ^ £• This implies Il^':,-/ll2 = \\PjU-h)h <e for ah j G Z. By the Riesz basis assumption on (j), we further have oo oo AY,\{fA,,k)? < WPjfWl < BY,\{f,^j,k)\\ (4.14) k=l k=l 13 c.f. Lemma 2.7 in jl8j . Further, if supp/ G [0, i?] then R {f,4>j,k) = {f,$,,k)= / /(A)p?(A)0(2-U)d^(A) Jo where Note that the functions {pl^\ A; € N} form an orthonormal basis of i^([0, 2-^], u;|[o,2j]) =• -^j- Suppose now that j is sufficiently large, i.e. 2^ > R. Then {f,<t>j,k) = (/?(2-^-),p1''^)x,. Thus by Parseval's equation for Xj, oo f-R EK/'-^^^)!' = l|/?(2"'-)|li. = / |/(A)n0(2-^A)pda;(A). k=i ' As is assumed to be continuous in 0, the functions A \4>{2^^X)\ converge to the constant \4>iG)\ > uniformly on [0,R] as j oo. Hence e > lim ||P,-/||2 > VA\m\\\fh > ^\m\i\\h\\2-e). As e is arbitrarily small, this shows that /i = and hence axiom (3) is satisfied. Vice versa, axiom (3) implies that lim Pjf = f for all / G L^{H). If / is compactly supported, then the same calculation as above shows that lim ||P,/||2 < VB\m\\\f\\2 j->oo which enforces ^(0) 7^ 0. Finally, the <j)j^k , k £'N are orthonormal if and only if it is possible to choose A = B = 1 in (|4.14|) . In this case we obtain (for / as just before) lim ||P,-/||2 = |0(O)|||/||2 = 1^(0)1 II/II2. Thus (3) is satisfied exactly if 10(0)1 = 1. □ Let us now write the two scale relation 1)4. 9() in a slightly different form, namely 0(2A) = G(A)0(A) (4.15) with G(A):^ sin(27rA) The filter function G is obviously 2-periodic and even. Whenever its restriction to [0, 1] is contained in -L^[0, 1] it can be expanded as a cosine series, 00 G(A) = V^y~^ ffn cos(n7rA). n=0 As for a classical multiresolution analysis one proves the following. 14 Lemma 4.10. Suppose that (j) G L'^{H) is an orthonormal scaling function of a radial MRA. Then the associated filter function G satisfies |G(A)|2 + |G(A + 1)|2 = 1 a.e.. (4.16) // in addition (j) £ (H) then ^4-16^ holds pointwise and G(0) = 1, G{1) = which implies oo oo V2^5n = 1, V2^(-ir5n = 0. n=0 n=0 Proof. In view of Corollary 14.31 we have oo oo 1= Yl l^(A + 2n)|2= |G(A/2 + n)|2|0(A/2 + n)|2 n=— OO n=—co oo oo = |G(A/2)|2^ |^(A/2 + 2n)|2 + |G(A/2 + l)|2 ^ |0(A/2 + 2n + 1)|2 n=— oo ?i=— oo = |G(A/2)p + |G(A/2 + l)p almost everywhere. If G (H) then (p is continuous and 0(0) 7^ by Theorem 14.91 Hence G(0) = 1 by (|4.15j) and G(l) = is an immediate consequence of (j4.16|) . □ 5 Orthogonal radial wavelets In this section we construct wavelets for a given radial MRA {Vj}ji^z in with orthonor- mal scaling function and filter function G. As usual, the wavelet space Wj is defined as the orthogonal complement of Vj in Vj+i, Thus LP'^H) decomposes as an orthogonal Hilbert sum 00 L\H) = W,. j=-oo Recall the definition of 5 in Section 4 and the characterization of Vq = K/> according to Lemma 14.51 feVo^ /(A) = ^ 0(A) with pes. (5.1) Define Sq := {a £ S : a(A + 1) = — a(A) for almost all A} which is a closed subspace of S w.r.t. || -111,2 [g^ij. Then W^i = Vq V-i is characterized as follows: 15 Proposition 5.1. (i) Let f E L^{H). Then f E W_i ^ /(A) = G{X + 1) ^(A) for some a G So- A (ii) The mapping Sq W-i ^ fa with UX) := (2vr)i/4^G(A + l)0(A) is an isometric isomorphism. Proof. For (i), notice first that the Hankel transform is a unitary isomorphism of L'^{H), so W-i = VqQV-i. Rescaling of by the factor 2 and relation (I^TT^ imply that h £ L'^{H) is contained in V-i if and only if there exists some f3 G S such that MA)=M?,2A) = Mg„)J(„. Thus P G S corresponds to / S W^i according to ()5.1|) if and only if /?(2A) ^^^^ ^|^(^)|2^^(^) ^ for all /3 G 5. A ' ' A Up to a constant factor, the integral on the left equals /?(2A)G(A) /5(A) |(/)(A)p dA = V / /3(2A) /3(A + 2n)G(A + 2n) |0(A + 2n)|2dA + V /" /3(2A) /3(A + 2n + 1)G(A + 2n + 1) |^(A + 2n + l)|2dA n&L' -1 /3(2A)(/5(A) G(A) + /3(A + 1) G(A + \))dX, where we used the periodicity and symmetry properties of I3.,I3^G as well as Corollary 14.31 Since /? G S" is arbitrary, we conclude that the vectors (/3(A), /3(A+1))"^ and (G(A), G(A+1))"^ must be orthogonal in for almost all A. This means that for some function a : [0, 1] C. Thanks to the boundedness of G, a belongs to -L^[0, 1]. Since [3 and G are 2-periodic on R, (3 is odd and G is even, an extension of a to M must be 2-periodic and odd; hence a G 5 with /3(A) = a(A)G(A + 1). Using the 2-periodicity of (3 and G in the second component of 1)5. 2|1 . we further deduce that /3(A) = — a(A + 1)G(A + 1) and therefore a(A + l) = — a(A) apart from the zero-set of G(A-|-1). Thus / is of the claimed form. Conversely, if a G 5*0, then /3(A) := a(A)G(A + 1) G S*, and (|5.2|) is satisfied. For the proof of (ii), we calculate Wf^Wl = wig =V2^jJ|"M^|G(A + l)|2|^(A)pd^(A) |a(A)|2|G(A + l)|2|^(A)|2dA = j' |a(A)|2|G(A + l)|2dA, 16 where we used that a and G are 2-periodic and </> is an orthonormal scahng function. By assumption on q, we have a(A — 1) = a(A + 1) = — a(A). Thus by Lemma I4.1U1 1 ,\ 1-1 |a(A)p|G(A + l)|2dA = / |a(A-l)|2|G(A)|2dA + / \a{X)\'^\G{X + l)\'^dX 1 1 |a(A)pdA = \\a\\l. This proves (ii). □ It is now easy to obtain an orthonormal basis of W-i. Recall that the Sfc, A; G N form an orthonormal basis of S. Moreover, the S2k~i, k £ N are an orthonormal basis of 5o. Thus by the previous result, the functions constitute an orthonormal basis of W-i. Define ip G L^{H) by V'(2A) = G(A + 1)(/>(A). (5.3) Then in view of (|4.3() . To obtain an orthonormal basis of Wq, we just have to rescale. Extending the notation r('^) := to r G N/2 and using the relation DaT^ = TaxDa, we obtain that an orthonormal basis of Wq is given by the functions V'. := Dy,.h = M^T^i^^-^m^, k G N. We call ij) a (basic) wavelet for the radial mutiresolution iVj). Definition 5.2. For j G Z and A; G N, define the "radial" wavelets V.,. ,(.) := D,-Mr) = 8^V2^!^|zir((2'=-i)/2)^(2i,). We have proven: Theorem 5.3. (i) For each j G Z, the set {ipj^k : A; G N} constitutes an orthonormal basis of Wj . (ii) The set {ipj,k : j G /c G N} is an orthonormal wavelet basis of L'^{H). Corollary 5.4. The functions ^j^k{x) := ^Pj,k{\x\), xeR^j ez,ken form an orthonormal basis for the closed subspace L^^^(M^) := {/ G L^(M'^) : / radial} of radial functions in L?' 17 6 Construction of radial scaling functions and wavelets Yet, we do not have handsome criteria in order to decide whether a given function (j) G Lp'{H) is a radial scaUng function, i.e. a scaUng function for a radial MRA. The analogy of our constructions to those on the group (M, +), however, leads to the following close relationship: Theorem 6.1. Suppose (p^ is a classical scaling function on M which is even and such that its (classical) Fourier transfrom is continuous in and satisfies J-{(j)u) G L^{H). Define 4> € -^^(^) ''^^'^ ^^-s Hankel transform, </.(A) := V2^.^(0k)(^A). (6.1) Then (j) is a radial scaling function. Conversely, if (j) is a scaling function for a radial MRA such that <j) is continuous in 0, then (j) £ L^(M) and the function cp^ defined by id.l]) (where (p is extended to an even function on M) is a classical scaling function on M. Moreover, (p is an orthonormal radial scaling function if and only if (p^ is an orthonormal classical scaling function. Proof. Let us start with the first assertion. As (p^ is a classical scaling function, we have by eq. (5.3.2) in fcez with suitable constants {) < A < B < oo. Moreover, c^r is orthonormal if and only if A = B = 1. Since (p^ is assumed to be even, definition ()6.1() is compatible with the even extension of <p. By Proposition 14.21 the set {Mj^T^^^(p, k E N} forms a Riesz basis for Vq = span{i?(^} which is an orthonormal basis if and only if (p^ is orthonormal; c.f. Corollary 14.31 Moreover, by eq. (5.3.18) of there exists a 27r-periodic function rriQ £ L^([— 7r,7r]) such that J^{<pM.){(,) = "T'o(C/2).^((/'r)(^/2), and mo is necessarily even in our case. Hence, with 7(A) := mo(A7r) sin(27rA) which clearly is contained in S, we have sin(27rA)(/)(2A) = 'y{X)(p{X). This is exactly the radial two-scale equation ()4.9|) . As (p is continuous in 0, the condition 0(0) 7^ of Theorem 14.91 is automatically satisfied (see e.g. Remark 3 on p. 144 in 6 ), and thus we finally obtain that (p is a radial scaling function. For the converse part notice first that continuity of i?!> in already implies that p £ (M) . We further proceed similar as before, using Proposition 5.3.1 and Proposition 5.3.2 in |S] and the corresponding results of the present paper. □ This theorem supplies a variety of radial scaling functions since there are many classical scaling functions on R which satisfy the assumptions of the theorem. However, as to or- thonormal radial scaling functions with compact support, a famous theorem of Daubechies implies the following negative result. Corollary 6.2. There do not exist any real-valued orthonormal radial scaling functions with compact support. Proof. The proof of Theorem 8.1.4 in shows that an even, real-valued and compactly supported scaling function is necessarily the Haar function X[-i/2,i/2]) the characteristic function of the interval [—1/2,1/2]. However, its Fourier transform J^{x[-i/2,i/2]){0 = ./ ig not contained in L^{H). □ 18 If (f) correponds to an even classical scaling function (f)^ according to Theorem 16 .11 then the hypergroup translates 0o,A: = ^kTk-w'P may be expressed accordingt to the formula Mx) = ^^-i^ (<^k(^ -k)- <Am(^ + k)) . (6.2) In fact, by the Plancherel theorem for the Hankel transform and eq. (|2.5|1 . ) -1 poo j{\r)j{\s)4>{\)du{\) = -j= / ^(A) sin(sA) sin(rA)dA V27r • rs J-oo for all r, s G M+. Here (p is as usual extended to an even function on M. Using relation H6.1|) as well as basic trigonometric identities and the Plancherel theorem for the classical Fourier transform, we can write Tr4){s) = — j= / (^(A)(cos A(r — s) — COS A(r + s))(iA 2\/27f • rs 1 /, (r — s\ , (r + s 27r • rs V V vr / V tt This implies (|6.2() . As an example, we consider the radial analogue of the Shannon wavelets. We define the scaling function via its Hankel transform, 2i\\ (W \ /2" sin(a;) -xcos(x) 0(A) = X[o,i](A), 4>{x) = ^3 . Constructing the associated basic wavelet according to formula (|5.l-i|) yields (after a short calculation) ; /2 sin(2x) — sin(x) — 2x cos(22;) + X cos(x) V'(A) = X[i,2](A), ^{x) = ^3 . The translates of the scaling function and the wavelet turn out to be , , 1 /sin(x - /cvr) sin(x + /cvr) \ , , , ^ 1 ( ( 2A:-1 \ / 2A:-1 \\ V'0,fc(2;) = ,^ ,1/, p X — TT -p\x^ TT , A; G N, with ^{x) sin(2x) — sin(a;) X 7 Algorithms For the use of our radial multiresolution in applications we need to formulate decomposition and reconstruction algorithms. The first step in such an algorithm consists of projecting the function / into a scale space Vj for some suitable j. We obtain a representation oo k=l 19 So from now on we assume that we have given a function / E V^- in terms of its coefficients (i) cj^ . The decomposition algorithm consists of decomposing / into Vj-i and Wj-i, i.e. of (j—i) (7 — 1) calculating the coefficients and dj^ in the representation 00 k=l k=l (Such a representation exists, since by construction {^j-i^fc, '>Pj-i,k '■ A; G N} is also a basis of Vj.) A reconstruction algorithm determines the coefficients cj^^ when / is given in terms of cjf "^^ and djf "^^ keN. We still assume that (j) is an orthonormal scaling function (and ip is hence an orthonormal wavelet). Let ■= {4>i,k,4>0,l) = {4'j,k,4'j~i,l), := (01,^,-00,/) = {4'j,k,il^j-i,l)- By using Hilbert space techniques - in particular Parseval's equation - we obtain analogously as in standard wavelet theory the decomposition formulae 00 00 ^£ ~ Z^^fc ' ~ 2-^^k 'i ' fc=i fc=i and the reconstruction formula 00 00 £=1 i=\ It turns out that the coefficients qf"^ and r^*"^ are determined in terms of the numbers gn in the cosine expansion of G, i.e. the coefficients in 00 G(X) = \/2 gn cos(n7r A) . n=0 Theorem 7.1. For £,k eN it holds (k) 1i gk-2e - 9 2£+k for 21 < k, 250 - 9M for 2£ = k, , g2e-k - g2i+k for 2£ > k, {-lf~^{gk^2i+i - 9k+2i-i) for 21 -Kk, 2go - 9U-2 for 2i-l = k, {-l)''^\g2e^i-k - 92£-i+k) for 2^ - 1 > A;. Proof. Using the Plancherel theorem, relation 1)4. 15(1 and Corollary 1)4. 3|) we obtain ={4'i,k,4'o,e) = {(i>i,k,(i>o/) "1/2 / n,(\l'}\n„(\l9'\n(\l9\\^(\l9'\\'^, Pk{\/2)p,{\/2)G{\/2)\(t){\/2)\'du:{\) ■V2 Sk{\)si{2\)G{X) \^{\ + 2n)\^dX •^^ n=-oo 00 „i -Ay^g^ / sin(A;A7r) sin(2£A7r) cos(n7rA)(iA. (7.1) n=0 20 An easy calculation using trigonometric identities shows sin(/i:A7r) sin(tA7r) cos(n7rA)dA t-k\ 3(25o,t-fc ^n,t+k) for n > 0, (5o,t+fc) for n = 0. (fe) (7.2) Setting t = 21 and inserting into (|7.1() yields the assertion for . We proceed similarly for rf =(</'i,fc,V'o/) = ^ / Sfc(A)s2£-i(A)G(A + l)dA JO sin(A;7rA) sin((2£ - l)7rA) cos(n7r(A + l))dA oo „i 4V5„(-1)" / sin(/fc7rA)sin((2£- l)7rA)cos(n7rA)(iA. oo n=0 oo Setting t = 2£ — 1 in ()7.2() and inserting into the last expression gives the result for r'l^ . □ Let us consider the case where only finitely many coefficients are different from zero. Although this is not possible for real-valued orthonormal scaling functions this assumption makes it easier to compare the radial wavelet algorithm with the classical one. Of course, in applications one can only handle finitely many coefficients anyway. So let us assume suppg C [0,A^], i.e. = for k ^ {0, ...,A^}. Elementary considerations show the following. Leaving k fixed yields (k) U I for£^[^,^] for £ ^ [1, ^] for £ ^ [^^, for ^ ^ [1, ^l±k±l] ifk>N, iik<N, i{k> N + 1, iik<N + l. If £ is fixed then (fc) U (k) u Sk) I With for /fc ^ [2£ -N,2i + N] iov k ^[l,N + 21] for /c ^ [2^ - 1 - iV, 2i - for A; ^ [1,2^ - 1 + iV] g\^ for 1 < \k\ 2gQ for k = 0, otherwise. if 2i > N, if 2i < N, l + N] if 2^ - 1 > iV, if 2£ - 1 < TV. <N-1, it holds G{X) = ■^'Ylik=-N^k(^^^^ ■ Because of the conditions on it, G is also the filter function for an ordinary multiresolution analysis on M with coefficients /i^. Now, if 2^ > resp. 2^ — 1 > A^ then it is easy to see that =hk for A; = -A^, . . . , N, ^fl-i+k) for fc = -A^,. 21 Similarly, ii k > N then ) = h^e, = {-lf-+'hk+i-2e. Hence, for 2i > N + 1 the decomposition formulae become oo oo k=l k=l and for A; > iV + 1 the reconstruction formula is oo oo £=1 £=1 These formulae are well-known. Indeed, they are the decomposition and reconstruction formulae of the classical discrete wavelet transform. So our approach leads to the classical algorithm if we are far enough away from the origin. If we are close to the origin we have derived an algorithm to handle the boundary point 0. References [1] J. P. Antoine, P. Vandergheynst, Wavelets on the n-sphere and related manifolds, J. Math. Phys. 39 (1998) no. 8, 3987-4008. [2] L. Baggett, H. Medina, K. Merrill, Generalized multi-resolution analyses and a construction procedure for all wavelet sets in R"., J. Fourier Anal. Appl. 5 (1999), 563-573. [3] D. Bernier, K.F. Taylor, Wavelets from square-integrable representations, SIAM J. Math. Anal. 27 (1996), 594-608. [4] W.R. Bloom, H. Heyer, Harmonic Analysis of Probability Measures on Hypergroups, de Gruyter, Berlin - New York, 1995. [5] M. Conrad, J. Prestin, Multiresolution on the Sphere, in: Tutorials on Multiresolu- tion in Geometric Modelling, Springer- Verlag, Berlin 2002. [6] I. Daubechies, Ten Lectures on Wavelets, SIAM, 1992. [7] J. Epperson, M. Frazier, An almost orthogonal radial wavelet expansion for radial distributions, J. Fourier Anal. Appl. 1 (1995), 311-353. [8] A. Erdelyi, W. Magnus, F. Oberhettinger, F.G. Tricomi, Higher Transcendental Functions, vol. U New York, Mc Graw Hill, 1953. [9] W. Freeden, T. Gcrvens, M. Schreiner, Constructive Approximation on the Sphere: with Applictions to Geomathematics, Clarendon Press, Oxford 1998. [10] J. Hinz, Hypergroup Actions and Wavelets, in: Infinite Dimensional Harmonic Anal- ysis (Japanese-German Symposium Kyoto 1999); eds. H. Heyer et al; Grabner- Verlag 2000, pp. 167-176. 22 [11] R.I. Jewett, Spaces with an abstract convolution of measures, Adv. in Math. 18 (1975), 1-101. [12] J.F.C. Kingman, Random walks with spherical symmetry, Acta Math. 109 (1965), 11-53. [13] W.R. Madych, Some elementary properties of multiresolution analyses of L'^{R^). In: Wavelets. A tutorial in theory and applications. Edited by Charles K. Chui. Wavelet Anal. Appl. 2, Academic Press, Boston, MA 1992, 259-294. [14] S.G. Mallat, Multiresolution approximation and wavelet orthonormal bases ofL'^{R.), Trans. Amcr. Math. Soc. 315 (1989), 69-87. [15] H. Rauhut, Wavelet transforms associated to group representations and functions invariant under symmetry groups, Preprint, 2003. [16] K. Trimeche, Generalized wavelets and hypergroups, Gordon and Breach, 1997. [17] G.N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, 1966. [18] P. Wojtaszczyk, A Mathematical Introduction to Wavelets, Cambridge University Press, 1997. Holger Rauhut Zentrum Mathematik Technische Universitat Miinchen D-80290 Miinchen Germany rauhut @ma.tum . de Margit Rosier Mathematisches Institut Universitat Gottingen Bunscnstr. 3-5 D-37073 Gottingen Germany roesler@uni-math.gwdg.de 23